Optical properties of small polarons from dynamical mean-field theory. 
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The optical properties of polarons are studied in the framework of the Holstein model by applying 
the dynamical mean-field theory. This approach allows to enlighten important quantitative and 
qualitative deviations from the limiting treatments of small polaron theory, that should be considered 
when interpreting experimental data. In the antiadiabatic regime, accounting on the same footing 
for a finite phonon frequency and a finite electron bandwidth allows to address the evolution of the 
optical absorption away from the well-understood molecular limit. It is shown that the width of 
the multiphonon peaks in the optical spectra depends on the temperature and on the frequency in 
a way that contradicts the commonly accepted results, most notably in the strong coupling case. 
In the adiabatic regime, on the other hand, the present method allows to identify a wide range 
of parameters of experimental interest, where the electron bandwidth is comparable or larger than 
the broadening of the Franck-Condon line, leading to a strong modification of both the position 
and the shape of the polaronic absorption. An analytical expression is derived in the limit of 
vanishing broadening, which improves over the existing formulas and whose validity extends to any 
finite-dimensional lattice. In the same adiabatic regime, at intermediate values of the interaction 
strength, the optical absorption exhibits a characteristic reentrant behavior, with the emergence 
of sharp features upon increasing the temperature — polaron interband transitions — which are 
peculiar of the polaron crossover, and for which analytical expressions are provided. 



I. INTRODUCTION 

The motion of electrons in solids is often coupled to 
the lattice degrees of freedom. There are broad classes of 
compounds where the electron-lattice interaction is such 
that the carriers can form small polarons, i.e. they are 
accompanied by a local lattice deformation that strongly 
modifies their physical properties, and can lead to the 
self-trapping phenomenon. This manifests in several ex- 
perimentally accessible quantities 1 . First, the d.c. con- 
ductivity is thermally activated in a broad temperature 
range, with a gap related to the energy of the electron- 
lattice bound stated. This is the energy barrier that 
the particle has to overcome to hop from site to site, 
which is proportional to the coupling strength. Second, 
most of the single-particle spectral weight — a quan- 
tity which is experimentally accessible through tunnel- 
ing or photoemission spectroscopy — moves to high en- 
ergies, corresponding to multi-phonon excitations in the 
polaron clou d 4 i 5 i 6 i 7 i 8 . Related to this, a broad contribu- 
tion emerges in the optical conductivity, typically in the 
mid-infrared region, which can be ascribed to transitions 
inside the polaron potential-wel&2ii2iiiii2ii£. 

Besides the strength, the range and specific mechanism 
of electron-phonon coupling, the phenomenon of polaron 
formation depends on the properties of the host lattice, 
such as the dimensionality, the width of the conduction 
band and the frequency of phonon vibrations. In ex- 
tremely narrow band systems, for example, the phonon 
energy can be of the order or even larger than the elec- 



tronic bandwidth. This leads to the concept of antia- 
diabatic quasiparticles, that can move through the lat- 
tice being accompanied by a very fast phonon cloud, 
which is at the very basis of the standard small polaron 
treatments 2 '^ 4 . In this framework, the small polaron op- 
tical absorption is strongly reminiscent of the behavior of 
electrons in a gas of independent molecules, and consists 
of a series of narrow peaks at multiples of the phonon fre- 
quency, whose distribution is determined by the electron- 
phonon coupling strength. 

However, narrow band materials exist where the 
phonon frequency is appreciably smaller than the elec- 
tronic bandwidth, which therefore cannot be described 
within the antiadiabatic approximation. In this case, 
one has to face the effect of the relatively large trans- 
fer integrals between molecules, which eventually wash 
out the discrete nature of the excitation spectra. Con- 
trary to the antiadiabatic situation, where polaron for- 
mation occurs through a smooth crossover, in the adia- 
batic case a sharp transition separates the weak-coupling 
regime from the polaronic regime in dimensions greater 
than li^iSiii. In the weak coupling limit, the optical ab- 
sorption at low temperatures consists of an asymmetric 
band resulting from the excitation and absorption of a 
single phonon. In the strong coupling regime, the ab- 
sorption spectrum moves to higher energies, reflecting 
the localized nature of the polaron. In this regime, pro- 
vided that the free-electron bandwidth is small compared 
to the broadening of the Franck-Condon line, the optical 
absorption acquires a typical gaussian lineshapeAiSi In 
the opposite situation, i.e. when the electron dispersion 
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dominates over the phonon-induced broadening, a differ- 
ent absorption mechanism sets in, corresponding to the 
photoionization of the polaron towards the free-electron 
continuum^2iI£ To our knowledge, there is at present no 
satisfactory theory of the optical spectra in the regime 
where both mechanisms coexist, which is often the case in 
real experimental systems. Another important unsolved 
question, that cannot be addressed by standard methods, 
is how the optical absorption evolves between the weak 
and the strong coupling limit, and recent calculationsiSiiS 
have shown that the spectra in the intermediate coupling 
regime exhibit specific features that are characteristic of 
the polaron crossover region. 

The simplest model which describes the rich phe- 
nomenology depicted above is the Holstein model^, where 
tight-binding electrons are coupled locally to disper- 
sionless bosons. The Holstein polaron problem can be 
solved^ in a non-perturbative framework using the dy- 
namical mean field theory*^ (DMFT). While this method 
neglects the spatial dependence of the electron self en- 
ergy, it treats exactly the local dynamics, which makes it 
particularly well suited to describe small polaron physics. 
It goes beyond the standard analytical approaches, as it 
can deal on the same footing with the low energy effective 
quasiparticles and the high energy incoherent features, 
which is crucial to the understanding of the polaron for- 
mation process^. The effects of a nonvanishing electron 
bandwidth and a nonvanishing phonon frequency are nat- 
urally included, which provides a unified description of 
the optical properties of small polarons, with no restric- 
tions on the regimes of parameters. Compared to fully 
numerical methodsiiiSiSl, which are also able to tackle 
the intermediate coupling regimes, the present approach 
is advantageous for the calculation of the optical conduc- 
tivity, since it is free of finite size effects, and gives direct 
access to the electronic excitation spectrum in real fre- 
quency. In contrast, finite cluster diagonalization studies 
suffer from the discretization of the Hilbert space, which 
can be quite severe in the polaronic regime where many 
phonon states are needed, while quantum Monte Carlo 
treatments usually work in imaginary time, and rely on 
analytical continuation algorithms for the extraction of 
spectral and optical properties. Note however that the 
present DMFT results, which are based on an exact con- 
tinued fraction expansion in the low density limit £ can- 
not be directly generalized to finite densities. 

The aim of this work is to take advantage of the DMFT 
to address the regimes of parameters not covered by the 
standard formulas of small polaron theory. Applying 
this non-perturbative method, we can critically exam- 
ine the ranges of validity of the usual limiting approxi- 
mations, and point out the quantitative and qualitative 
discrepancies arising in several regions of the parame- 
ters space. Special emphasis is given to the following 
points, which are not accessible by the usual methods 
available in the literature: i) the evolution of the multi- 
peaked spectra in the antiadiabatic regime, at finite val- 
ues of the free-electron bandwidth; ii) the effects of a fi- 



nite electron dispersion on the usual Franck-Condon line- 
shapes, in the adiabatic polaronic regime; iii) the pecu- 
liar features arising in the region of the adiabatic polaron 
crossover. Concerning the last point, the results of our 
previous work Ref^S, where such features were first re- 
ported, are here extended to much lower temperatures, 
which is made possible by a new adaptive method in- 
troduced to deal with the fragmented excitation spectra 
characteristic of the polaron problem^ This new proce- 
dure allows to identify a reentrant behavior of the optical 
properties in the crossover regime, where increasing the 
temperature switches from a weak-coupling like absorp- 
tion to a typical polaronic lineshape. Furthermore, a 
previously published formula for the T = absorption in 
the polaron crossover regime is corrected here, and ana- 
lytical expressions are derived to describe the lineshape 
of the polaron interband transitions reported in Ref^. 

The paper is organized as follows: In Section II we in- 
troduce the DMFT formalism and the Kubo formula for 
the optical conductivity, and discuss the details of the 
calculations. Sections III and IV are devoted to the anal- 
ysis of the absorption spectra respectively in the antiadi- 
abatic and in the adiabatic regimes. The specific features 
arising at intermediate values of the coupling strength in 
the adiabatic polaron crossover region are examined in 
Section V. The main results are summarized in Section 
VI. 



II. MODEL AND FORMALISM 

We study the Holstein hamiltonian, where tight- 
binding electrons (c,- lCT , c] ) with hopping amplitude t are 
coupled locally to Einstein bosons (a,, a\) with energy lvq: 

H = uj ajai c+ a c l ^(a+ + ai) 

i i,cr 

- *£0*W + h.c.). (i) 

The single polaron problem can be solved in the frame- 
work of the DMFT 2 °, which yields an analytical expres- 
sion for the local self-energy E(j/) in the form of a con- 
tinued fraction expansion^. The latter must be iterated 
in an appropriate self-consistent scheme, where the free- 
electron dispersion defined by the tight-binding term in 
Eq. Q enters only through the corresponding density 
of states (DOS) N(e). From the knowledge of the self- 
energy, one can define the single particle spectral function 

v) = --Im ^j- (2) 

that carries information on the spectrum of excited 
states, and its momentum integral, the spectral density 
N*(y). The latter can be evaluated by introducing the 
Hilbert transform of the DOS, H(z) = J deN(e)/(z - e) 
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N*(v) = -(1/tt) Im H[v - E(i/)]. 



(3) 



The conductivity at finite frequency is related to the 
current-current correlation function through the appro- 
priate Kubo formula. In DMFT, due to the absence of 
vertex corrections^ * 2 ^, such two-particle response func- 
tion can be expressed as a functional of the fully inter- 
acting single particle spectral function p(e, v) (or, equiv- 
alently, of the local self-energy £). This can be related 
to the single polaron solution of Refi by performing 
an expansion in the inverse fugacitjii^ z~ x = exp(ju/T). 

— > as the chemical potential [i — ► -co in the limit 
of vanishing density at any given temperature). The opti- 
cal conductivity turns out to be proportional to the car- 
rier concentration x, and can be expressed in compact 
form as 



(1 



_ V(u,(3) 



(4) 



where the constant C = e 2 a 2 /Tiv carries the appropriate 
dimensions of conductivity (a being the lattice spacing, 
v the volume of the unit cell) and (3 — 1/T is the inverse 
temperature. We have defined 



V(fi>,0) = J deN(e)4,(e) x (5) 
x J dve-^ v - E ^p{e,v)p{e,v + uj) 
Af(j3) = f dve-^- Eo) N*{v) (6) 



where E is the polaron ground state energy at T = 0. 
The exponential factors account for the thermal occu- 
pation of the electronic levels (the Fermi temperature 
vanishes in the low density limit, and the particles obey 
Maxwell-Boltzmann statistics). The normalization fac- 
tor Af represents the partition function of the interacting 
carrier. Unless differently specified, the density of states 
N(e) of the unperturbed lattice, that enters in the def- 
inition of the correlation function T>, is assumed semi- 
elliptical of half bandwidth D 



N(e) 



tD 2 



Vd 2 



(J) 



corresponding to a Bethe lattice in the limit of infinite 
connectivity. This choice reproduces the low energy be- 
havior of a three-dimensional lattice, while leading to 
tractable analytical expressions. The function <j>{e) is the 
corresponding current vertex 



(D 2 



! )/3 
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that can be derived by enforcing a a sum rule for the to- 
tal spectral weightj2&24. Compared to the more rigorous 
procedure of Refi 2 ^, this has the advantage of leading to 
the correct threshold behavior for the optical conductiv- 
ity expected in a three-dimensional system in the weak 



coupling limit ^ as well as in the adiabatic photoioniza- 
tion limit as will be shown below. In the usual 
strong coupling regime, 9,10 on the other hand, the re- 
sults become independent on the detailed shape of both 
N and cj>, and the choice of the vertex function is not in- 
fluent. Let us mention that low-dimensional systems can 
also be studied with the present method, by replacing 
N(e) and 4>(e) of Eqs. Q and iJSJ with the appropriate 
DOS and vertex function. In such cases, the neglect of 
vertex corrections in the current-current correlation func- 
tion that is implicit in the DMFT formulation could in 
principle constitute an important limitation. In practice, 
however, the limiting cases of weak and strong coupling 
are correctly reproduced by the DMFT, which gives some 
confidence on the results obtained in more general situa- 
tions (see e.g. the discussion at the beginning of Section 
IV A). 

With the present choice of Eqs. Q) and JSJ, the Hilbert 
transform JC(z) of the product N(e)<p(e) can be expressed 
in closed form, which allows the e-integral in Eq. JSJ to 
be performed explicitely: 
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with 



dve- f3iu - Eo) B[oj + v -E(u + u),v- £(»] (9) 



B{z 1 ,z 2 ) = --^Re 



JC( Zl ) - JC(z 2 ) K{zt)-K{z 2 ) 



z 2 - Zl 
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The remaining integral in v appearing in Eq. (JHJ has 
to be computed numerically, which presents two main 
difficulties. First of all, at temperatures T « wo, the 
function B cx JmS(v) decays exponentially at energies 
v below the ground state E . Nevertheless, such regions 
contribute to the final result due to the presence of the 
thermal factor e~P( u ~ E °\ As a rule of thumb, at low 
temperatures, the integration limits must be extended 
down to v — Eq ~ —2Ep to account for a proper num- 
ber of excited states^. At the bottom of the integration 
region, the number B is therefore multiplied by an expo- 
nentially large factor e 2E p/ T which amplifies the numer- 
ical errors and, for typical values of Ep/loq ~ 10, limits 
the attainable low temperature limit to T/loq ^0.1. The 
second difficulty is that, in the strong coupling regime, 
the function B is zero almost everywhere for v < Eq, be- 
ing concentrated in exponentially narrow peaks of width 
~ e - £ p/"o ^ corresponding to multi-phonon resonances at 
multiples of wo (see Fig. below). For example, taking 
the parameters of Fig. [31b, the ratio between the width of 
the peaks and their separation is 10 -3 . Special care must 
be taken to manage with this rapidly varying function. A 
uniform mesh discretization of the integral appearing in 
Eq. © was used in RefpiS, which was appropriate for the 
intermediate coupling/temperature regimes. To attain 
the strong coupling and low temperature regimes, we use 
here an adaptive non-uniform mesh optimized to account 
for the narrow peaks of the function e~^^~ E °^ ImT,(v) . 

Bearing these limitations in mind, we shall present in 
the following Sections the numerical results for coupling 
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strengths up to Ep/uiQ ~ 10 and temperatures down to 
T/u! ~ 0.01 — 0.1 (the lowest temperatures are reached 
in the weak/intermediate coupling regime). Analytical 
formulas will be derived to access the limits of strong 
coupling and vanishing temperatures. For simplicity, we 
shall drop the numerical prefactor xC in Eq. J2J and 
express all energies in units of the half-bandwidth D. 



III. OPTICAL CONDUCTIVITY IN THE 
ANTIADIABATIC REGIME 

The process of polaron formation at zero temperature 
is different depending on the value of the adiabaticity 
ratio 7 = loq/D, which measures the relative kinetic 
energies of phonons and band electrons (D is half the 
free-electron bandwidth, proportional to the hopping am- 
plitude t). In the antiadiabatic regime (7 ^> 1), the 
buildup of electron-lattice correlations occurs through 
a smooth crossover controlled by the coupling parame- 
ter a 2 = Ep/ojor^ which sets the average number of 
phonons in the polaron cloud (Ep = g 2 /oJo is defined as 
the polaron binding energy on an isolated molecule). In 
this case the behavior of the optical absorption at any 
value of the electron-phonon coupling strength can be 
deduced to a good approximation from the analysis of an 
isolated molecule, corresponding to the limit D — > 0^2*^ 
It consists of a series of narrow peaks at multiples of 
u>q, whose distribution is determined by the coupling pa- 
rameter a 2 , and evolves gradually through the polaron 
crossover located around a 2 ~ 1 (see also Fig. [21 in the 
following). At low a 2 , the spectral weight is mainly lo- 
cated in the Drude peak, plus a weaker absorption peak 
at lu = ujq. Upon increasing the coupling strength, multi- 
phonon scattering processes become important, and sev- 
eral peaks arise at multiples of luq. For a 2 3> 1, the 
distribution of the peak weights eventually tends to a 
gaussian centered at lu = 2Ep = 2a 2 u>o, with a variance 
~ V / 2E p ll> [cf. App. A and Eq. (12) below]. 

Thermal effects lead to some redistribution of spectral 
weight among the different peaks: In the weak coupling 
regime, the thermal excitation of phonons generates addi- 
tional absorption peaks at multiples oiu>Q. In the strong 
coupling regime, increasing the temperature broadens the 
distribution of peak weights: For T > Wq/2, the variance 
tends to \JAEpT and is completely determined by the 
thermal fluctuations of the phonons (cf. App. A); When 
the temperature is increased further, thermal dissocia- 
tion of the polaron at T ~ Ep eventually washes out the 
absorption maximum at finite frequency, resulting in a 
transfer of spectral weight towards u> = 0. 

The above picture is based on the limit of vanishing 
electron bandwidth. A small nonvanishing transfer in- 
tegral can be expected to give the individual peaks a 
finite width, without modifying the distribution of spec- 
tral weights^. In the weak coupling regime, this scenario 
is confirmed by the DMFT results illustrated in Figure 
for 7 = 4 and a 2 — 0.25. The absorption spectrum 



at T = 0.2wo consists of few peaks of width ~ AD, cen- 
tered at multiples of luq; Increasing the temperature leads 
to a moderate transfer of spectral weight to higher fre- 
quencies, while the individual peaks shrink as predicted 
by Holstein's approximation 2 . An interesting behavior is 
seen at extremely low temperatures, where T < D < wo 
(inset of fig l.a). In this limit, the optical absorption 
is dominated by transitions between "low momentum" 
states, i.e. states located near the bottom of the sub- 
bands. As a result, the peaks become asymmetric and 
the absorption thresholds move from lu = uluq — 2D to 
lu = uloq [details arc given in Appendix A, Eq. (|25|) and 

mi 

The optical conductivity in the strong coupling regime 
(7 = 4 and a 2 = 5) is shown in Fig. [Jib. In this case, 
although the overall picture deduced from the molecular 
limit is qualitatively recovered, the fine structure exhibits 
a behavior that contradicts the standard results available 
in the Iiterature2i£»i2ii4, according to which the absorp- 
tion spectra at nonzero D consist of a series of narrow 
peaks of equal width, proportional to the renormalized 
polaronic bandwidth [this follows, through Eq. Ipjl. from 
the assumption that the multiphonon subbands in the 
single particle spectral density all have the same width] . 
Instead, the multiphonon peaks in Fig. ^b are broader 
at high frequency than at low frequency, although this is 
hardly visible on the scale of the figure. More surpris- 
ingly, their width increases with temperature, as shown 
in the inset of Fig. [I]b, contrary to what could be 
expected 2 . 

This apparently anomalous behavior can be traced 
back to finite bandwidth effects on the excitation spec- 
trum, that go beyond Holstein's decoupling scheme, and 
have already been reported in Refii as well as in Ref. 3 , 
where they have been shown to strongly affect the abso- 
lute value of the d.c. conductivity. In fact, the broaden- 
ing of molecular levels induced by band overlap is not uni- 
form over the excitation spectrum, but rather increases 
as they move away from the groud state (see e.g. Figs. 
17 and 18 in Refi): while the width of the lowest peak, 
related to coherent tunneling between different molecular 
units, scales exponentially as exp(— a 2 ), the width of the 
higher order peaks is much larger, being determined by 
both incoherent hopping between neighboring molecules 
and coherent hopping in the presence of multiphonon ex- 
cited states: for such peaks, the DMFT results indicate 
a power law dependence in a 2 . Moreover, as the temper- 
ature increases, the different levels are mixed by phonon 
thermal fluctuations. As a result, the width of the sub- 
bands at low energy rapidly increases with temperature 
to attain the typical broadening of the high energy sub- 
bands, as seen in the inset of Fig. ^b (in the data at 
T = 0.2w , a fine structure due to the superposition of 
at least two contribtions of different widths is visible). 

Note that the discrete multi-peaked structure charac- 
teristic of the antiadiabatic regime remains at all tem- 
peratures, contrary to what is stated in Ref^S: when 
DCwo, the thermal broadening of the peaks is not suf- 
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FIG. 1: (Color online) Optical conductivity in the antiadiabatic regime, a) for 7 = uJq/D = 4 and a 2 — 0.25 and b) a 2 = 5. 
a) In the weak coupling regime, the individual absorption peaks shrink upon increasing the temperature above T ~ ujq/2, as 
predicted by Holstein's approximation. In the inset, the low-temperature evolution of the one-phonon absorption is shown (full 
red lines, from left to right, T/loq = 0.2,0.05,0.02). The black dotted line is the weak coupling result HIGH at T = 0. The 
absorption edge shifts from uo = u>o — 2D at T > D to ui = u>o at T <C D (see appendix A), b) In the strong coupling regime, the 
optical absorption consists of extremely narrow peaks, similar to the case of an isolated molecule (the peaks at T = 0.6u>o and 
T — 2u>o have been shifted laterally by a constant offset for clarity). Contrary to the weak coupling limit, here the individual 
peaks broaden with temperature (see inset). 



ficient to lead to a continuous absorption curve, even at 
T 3> Ep, unless a sizeable phonon dispersion is intro- 
duced in the model (this would give the individual peaks 
an intrinsic finite width, preventing the extreme band 
narrowing in the limit D — > 0). 

IV. OPTICAL CONDUCTIVITY IN THE 
ADIABATIC REGIME 

In the adiabatic limit (7 = 0), the carrier properties 
change drastically from weakly renormalized electrons to 
self-trapped polarons at a critical value of the coupling 
strength A = Ep/D, above which a bound state emerges 
below the bottom of the free-electron band. For example, 
in the case of the semi-elliptical density of states Eq. (7), 
adiabatic polaron formation takes place at A c = 0.843. 
Allowing for lattice quantum fluctuations changes this 
localization transition into a very sharp crossover (cf. 
Fig. [5J which separates weakly renormalized electrons 
for A ^ A c from polarons with very large effective mass 
at A > A c £±& 

To describe this behavior, methods such as the DMFT 
that do not rely on a "small" parameter are extremely 
valuable for the following reasons. On one hand, there is 
at present no unified analytical approach which allows to 
calculate the optical conductivity (even approximately) 
in the whole range of parameters, from the weak to the 
strong coupling limit. This is not surprising, since the 
physics is fundamentally different in the two regimes. 
This situation should be contrasted with the antiadia- 
batic case, where the basic qualitative features of the 
optical absorption can be inferred form the solution of 
the model on a single molecule (cf. beginning of Section 



III A). In addition, even within the adiabatic strong cou- 
pling regime, the polaronic lineshape changes depending 
on the ratio between the broadening s of electronic lev- 
els induced by phonon fluctuations, and the free electron 
bandwidth D. Theories exist for the limits s/D — > 
and s/D — > 00, but their validity is questionable in the 
intermediate range of experimental relevance. 

In the following, we shall first report the analytical ex- 
pressions recovered within the present DMFT formalism, 
in three limiting cases: weak coupling, strong coupling 
s/D — > 00 and s/D — > 0. The first two formulas turn out 
to be perfectly equivalent to the results available in the 
literature, showing that at least in such limits, the neglect 
of vertex corrections in the Kubo formula implicit in the 
DMFT approach has negligible influence on the results. 
The third formula, on the other hand, constitutes an im- 
provement over the expressions of Refs^2il£, even when 
applied to one-dimensional lattices. The DMFT results 
obtained in more general cases will be presented next, 
pointing out the inadeguacy of the standard descriptions 
in many cases of interest. 



A. Limiting cases 

In the weak coupling limit (A — > 0), the optical ab- 
sorption at T = consists of a broad band related to the 
excitation of a single-phonon, with an edge at to = ujq 
followed by a power law decay at higher frequencies, and 
eventually an upper edge at uo = u; + 2D [cf. Eq. if^BI) 
in Appendix A] 

a{u) = ^j^<f>(u - D - uo)N(u - D - uq). (10) 
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In a three dimensional system, the absorption edge be- 
haves as (w — lu ) 3 / 2 . As the temperature increases, the 
gap below the threshold is rapidly filled and the absorp- 
tion maximum is washed out above T ~ u>o/2. 

In the strong coupling regime (A — * oo), the photoex- 
citation of the electron is much faster than the lattice 
dynamics, which is virtually frozen during the absorp- 
tion process. Since the lattice energy cannot be relaxed, 
the dominant optical transition corresponds to the dif- 
ference in electronic energy between the initial and final 
states (Franck-Condon principle) which, in the Holstein 
model, equals twice the ground state energy 2Ep. The 
shape of the optical absorption will depend on the ratio 
between the width of the non-interacting band ~ D, and 
the variance s of the phonon field, which controls the 
broadening of electronic levels. The latter obeys^i 



= EpujQ cothwo/2T 



(11) 



(cf. Appendix A). It is determined by the quantum fluc- 
tuations of the phonons at low temperatures, and in- 
creases due to the thermal fluctuations as T > wq/2. 

When s ^> D 7 i.e. when the phonon induced broad- 
ening of electronic levels is much larger than the elec- 
tronic dispersion, the absorption by localized polarons 
takes the form of a skewed gaussian peak centered at 

Umax = 2Ep*^— 



a(cj) = exp 

a 4 



(lu - 2E P f 



As 2 



(12) 



Following Eq. the Franck-Condon line further 

broadens upon increasing the temperature above T ~ 
luq/2 and eventually moves towards lo = at temper- 
atures higher than the polaron binding energy, as the 
polaron thermally dissociates. Note that the above for- 
mula also describes the envelope of the discrete absorp- 
tion spectrum of polarons in the antiadiabatic regime, 
shown in the preceding Section (cf. Fig. Q]b). 

To recover the standard result for a three-dimensional 
cubic lattice^ where the total bandwidth is 2D = 12t, 
the above formula has to be multiplied by a prefactor 
2/9, which corrects for the larger value of the average 
square velocity in the DOS Eq. (7) used in the calcula- 
tions. Analogous prefactors should be included when any 
of the DMFT results obtained in this work are applied to 
finite-dimensional lattices, that can be straightforwardly 
obtained by evaluating the second moment of the corre- 
sponding non-interacting DOS. 

In the opposite limit s <C D, the lineshape is domi- 
nated by the electronic dispersion. The absorption is due 
to transitions from a polaronic state whose electronic en- 
ergy is ~ —2Ep to the continuum of free-electron states. 
We have [cf. Eq. in Appendix A] 

4E 2 1 — e~"/ T 
cf{uj)=tt^- ef,(w - 2E p )N(lu - 2E P ) (13) 



In three dimensions, the absorption vanishes as Aw 3 / 2 
at the edges, as in the weak coupling case of Eq. (10). 
Taking a semi-elliptical DOS as representative of a three- 
dimensional lattice, we also find that the absorption max- 
imum LUmax = 2Ep — D 2 /2Ep is shifted to lower frequen- 
cies compared to the usual estimate. This softening is en- 
tirely due to finite bandwidth effects. Note that Eq. i|13|) 
is valid at all temperatures below the polaron dissocia- 
tion temperature T ~ Ep. In particular, contrary to Eq. 
(12), nothing happens here at temperatures T ~ ujq/2, 
provided that the condition s -C D is not violated. 

Formulas similar to Eq. (|13|) have been derived in 
Refii£, in Ref. 13 for one-dimensional systems, and in 
Refiii for the high frequency absorption threshold of 
large polarons in the Frohlich model. As discussed in 
Appendix A, these derivations are strictly valid only in 
the limit Ep — > oo. On the other hand, the present for- 
mula accounts for the non-negligible dispersion of the 
initial state at finite values of Ep, which is reflected 
in the additional prefactor 4Ep/uj 2 , leading in general 
to a much more asymmetric lineshape. Remarkably, 
if the appropriate one-dimensional DOS is used, Eq. 
(|13|l describes much better the one-dimensional exact 
diagonalization data of Refji 3 " than their own Eq. (9) 
[also reported in Appendix A as Eq. The ab- 

sorption maximum in this case is located at u) max = 
2E P - E P {y/\ + 6(D/E P ) 2 - l)/2. 

The regions of validity of the analytical formulas pre- 
sented in this Section are summarized in Fig. [21 
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We see that in this case the absorption of photons is 
only possible in the interval 2Ep — D < lo < 2Ep + D. 
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FIG. 2: (Color online) Phase diagram illustrating the regimes 
of validity, at T — 0, of the different limiting formulas derived 
in the text. The bold and dotted lines indicate the location 
and spread of the polaron crossover as defined in Refi±&. The 
(red) shaded region delimits the antiadiabatic regime 7 > 1 
treated in Section III. In the hatched region at 7 < 1, the 
optical conductivity is not accessible by the usual methods, as 
pointed out in Section IV. The (gray) shaded region indicates 
the adiabatic intermediate regime described in Section V. 



B. DMFT results 

The DMFT results obtained in the adiabatic regime 
are illustrated in Fig. [3] First of all, our data show that 
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FIG. 3: (Color online) Optical conductivity in the adiabatic regime, a) 7 = u)q/D — 0.1 and A = 0.7 and b) A = 1.1. a) In the 
weak coupling regime, the gap below uj = ujo is rapidly filled upon increasing the temperature, and the one-phonon absorption 
threshold is washed out at T > ujq/2. The black dotted line is the analytical weak-coupling result illlH . At higher temperatures, 
the low-frequency absorption is reduced as some spectral weight is transferred to higher frequencies (see inset: from top to 
bottom, T/loo = 0.4, 0.6 and 1). b) In the polaronic regime, the absorption peak broadens at T > ujq/2 and then moves towards 
uj — as T > Ep. The inset compares the low temperature spectrum with the limiting lineshapes Eq. (12) (full line) and Eq. 
(13) (dashed line). 



the absorption shape l|l(J|) derived in the limit of A — * 
remains qualitatively valid at finite values of A, up to 
the polaron crossover. This is exemplified in the data 
reported in Fig. [31a. at 7 = 0.1 and A = 0.7, i.e. at 
a value of the coupling strength which lies a priori well 
beyond the range of validity of perturbation theory^ not 
far from the polaron crossover coupling strength A c (7 = 
0.1) ~ 0.93 (see Fig. Eland Section V below). Despite 
such relatively large A, the threshold at ujq and the overall 
behavior of the high frequency tail at low temperatures 
(full red line) are quite close to the predictions of the 
weak-coupling formula i|10fl (black, thin dotted line). The 
increasing importance of multi-phonon processes shows 
up mainly in the fine structure, with the appearance of 
alternating peaks and dips at multiples of ujq. 

Figure 0b shows the DMFT results at 7 = 0.1 and 
A = 1.1. This value of the coupling strength lies in the 
polaronic regime A > A c , in a region where the electronic 
dispersion and the phonon-induced broadening are com- 
parable (s/D = 0.33), so that none of Eqs. (12) and (13) 
is expected to hold. The comparison with such limiting 
formulas is illustrated in the inset of Fig. 0b. While 
the position of the absorption maximum seems to agree 
with the prediction of Eq. (13) (dashed line), the peak 
height is much reduced and the absorption edge is com- 
pletely washed out by phonon fluctuations. On the other 
hand, the temperature dependence qualitatively agrees 
with Eq. (12) (full line), i.e. the peak broadens as the 
temperature is raised above T ~ ojq/2, which can be un- 
derstood because the low frequency tails are dominated 
by the phonon fluctuations. We see that the DMFT spec- 
trum not only lies at lower frequencies compared to Eq. 
(12), but it is also much broader and asymmetric (the 
same trend was observed in Ref^S in the one-dimensional 



case). All these effects can be ultimately ascribed to the 
finite electron bandwidth D. 

Our results show that detectable deviations from Eq. 

(12) arise as soon as the non-interacting bandwidth is 
larger than the broadening s, a condition that is com- 
monly realized in real systems. For example, taking typ- 
ical values coo ^ 0.01 - 0.05eV and E p ~ 0.1 - 0.5eV 
yields a zero temperature broadening s ~ 0.03 — 0.16el^, 
in which case electron bandwidths of few tenths of eV 
are already sufficient to invalidate the standard gaussian 
lineshape Eq. (12). In the opposite limit, marked devi- 
ations from Eq. (13) appear already at small values of 
the ratio s/D, roughly s/D > 0.05: as soon as a small 
finite broadening is considered, exponential tails of width 
s arise that wash out the sharp absorption edge of Eq. 

(13) , while the peak height is rapidly reduced (this can 
be partly ascribed to a rigid increase of the renormalized 
band dispersion, as evidenced in Ref^). The intermedi- 
ate region where both Eq. (12) and (13) fail in describing 
the polaronic optical absorption is shown as a hatched 
area in Fig. |21 

As was mentioned above, the ability to describe the 
continuous evolution of the absorption shape between the 
limiting cases of Eq. (12) and Eq. (13) ultimately follows 
from the fact that finite bandwidth effects are properly 
included in the DMFT. On the other hand, accounting 
for a non-zero phonon frequency allows to address the 
multi-phonon fine structure of the optical spectra, which 
is seen to evolve gradually upon varying the adiabaticity 
ratio 7. In fact, the optical absorption of a polaron at 
finite values of 7 shares features with both the antiadia- 
batic and adiabatic limits, having a pronounced structure 
at low frequencies (possibly with multiple separate nar- 
row bands), followed by structureless tails at higher fre- 
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quencies (cf. Fig. |3|b). As 7 increases, the fine structure 
progressively extends to higher frequencies, and evolves 
into the discrete absorption pattern of Fig. ^b, typical 
of the antiadiabatic limit. 



V. POLARON CROSSOVER IN THE 
ADIABATIC REGIME 

The critical coupling strength for polaron formation 
was identified in Ref as the value at which oIEq / dg has 
maximum slope, which follows approximately A c (7) ~ 
0.84 + O.97 at small 7 (note that this derivative is re- 
lated, through the Hcllmann-Fcynman theorem, to the 
electron-lattice correlation function). The width of the 
crossover region, obtained by looking at the maximum 
slope of \d 2 E /dg 2 \, is roughly given by A A ~ 7/2 (cf. 
Fig. 3.b in refi&, reported in Fig. Inhere), and vanishes in 
the adiabatic limit where the polaron formation becomes 
a true localization transition. 

In this Section, we address the specific properties of 
the optical absorption in the adiabatic polaron crossover 
region. The present results extend the preliminary re- 
sults obtained in Ref. 18 allowing an inspection of the very 
low temperature regime, which turns out to be crucial to 
understand the evolution of the optical properties. We 
show that the adiabatic crossover region has two original 
signatures: The first is a coexistence of both weak and 
strong coupling characters in the optical conductivity, 
with a transfer of spectral weight between the two occur- 
ring at very low temperatures. The second is the emer- 
gence of narrow absorption features at low frequencies (of 
the order of, or even below the phonon frequency) with a 
non monotonic temperature dependence, corresponding 
to resonant transitions between long-lived states in the 
polaron excitation spectrum^* We finally provide a semi- 
analytical expression for the optical conductivity in the 
zero temperature limit, that corrects Eq. (7) of Ref»i£. 



A. Low temperature transfer of spectral weight 

A typical optical conductivity spectrum in the 
crossover region is shown in Fig. 0]a, for 7 = 0.1 and 
A = 0.9. The data at the lowest temperature (red full 
curve) are very similar to the A = 0.7 results of Fig. [21a, 
with a marked singlc-phonon edge at lu = loq as given by 
Eq. (|TU|l , followed by a multi-phonon structured high fre- 
quency tail. Such "weak-coupling" behavior is in agree- 
ment with the fact that A = 0.9 lies slightly below the 
critical value A c (7 = 0.1) ~ 0.93. Upon increasing the 
temperature (green dashed curve), the optical absorp- 
tion develops a broad maximum at high frequency typi- 
cal of the polaronic regime, comparable to what is seen 
at A = 1.1 in Fig. |3|b. Remarkably, this occurs at tem- 
peratures much below the phonon quantum to thermal 
crossover at T ~ wo/2, suggesting that such transfer of 
spectral weight is governed by an electronic energy scale. 



To understand the origin of this phenomenon, a more 
detailed analysis of the single-particle spectral function 
is needed. As can be seen in Eq. (JSJ, the calculation of 
the optical conductivity involves a convolution of the two 
quantities p(e, v + u>) and e~~P( v ~ E ° > p(e, v). Due to the 
exponential weighting factor, the contributions to p(e, v) 
at v > Eq are strongly suppressed at low temperatures, 
and only the states related to phonon emission processes 
at v < Eo will contribute. For such incoherent states, 
the spectral function is directly proportional to the scat- 
tering rate 7m£(z/) [cf. Eq. (JI5J in App. B]. There- 
fore, neglecting the e-dependence in Eq. (|SJ, the optical 
spectrum is roughly proportional to a convolution of the 
spectral density N*{v + ui) with the weighted scatter- 
ing rate e~^^ u ~ E °^ ImY,(y). Both of these quantities are 
illustrated in Fig. [SJ 

With the present choice of parameters, the spectral 
density at low temperature (upper panel of Fig. |SJa) con- 
sists of a single subband of width W < ujq (the polaron 
band, expanded in Fig. [SJc) disconnected from the elec- 
tron continuum at higher energy. Scattering processes 
are suppressed below Eq + ujq because we are consid- 
ering dispersionless optical phonons, so that the whole 
polaron band represents coherent long-lived states at 
low temperature. In contrast, higher energy states at 
v > Eq + ujq are mostly incoherent, being strongly scat- 
tered by phonons. The marked asymmetry of the po- 
laron band in Fig. [5Jc comes from the fact that the band 
dispersion flattens in proximity of the band top, where 
the states have a more localized character. This coex- 
istence of free-electron like states at low momenta and 
localized states at high momenta 17 is crucial here, and 
can be understood from the following argument: In the 
weak coupling limit, the spectral density develops a dip 
at v = Eq + luq corresponding to the threshold for one- 
phonon scattering. Below this value, the band dispersion 
flattens due to the hybridization with the (dispersionless) 
phonon states, leading to a very large density of states. 
When A increases, the dip eventually becomes a true gap 
and a polaron band of width W < loq emerges from the 
continuum, retaining a characteristic asymmetric shape 
with a maximum close to the top edge. 

The excitation spectrum hardly changes when going 
from T = 0.03^0 to T — 0.2cj , except for the appear- 
ance of few exponentially weak replicas of the polaron 
band at v < Eq (see Figs. |SJa andjSJc). On the contrary, 
the weighted scattering rate (lower panel of Fig. |5]a) un- 
dergoes a drastic change in the same temperature range. 
At T = 0.03wo, the latter consists of a pattern of equally 
spaced peaks at v = Eq — vicq, whose intensity decays as 
v moves away from the ground state energy. Retaining 
only the (largest) peak at Eq — ujq, we see that the convo- 
lution integral T> roughly follows the form of the spectral 
density N*(lo + luq — Eq), with a gap below uj = luq and 
a maximum around uj « Ep. The behavior changes at 
T = 0.2wo, where the weighted scattering rate acquires a 
gaussian distribution centered at an energy ~ Ep below 
the ground state, so that the maximum of T> moves to 
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FIG. 4: (Color online) Optical conductivity in the adiabatic regime 7 = 0.1, at an intermediate value of the coupling strength 
A = 0.9. a) the absorption is "weak-coupling like" at low temperatures, with a clear threshold at u) = u>o, and becomes 
"strong-coupling like" at higher temperatures, causing a non-monotonic evolution of the spectral weight, b) a detailed view of 
the temperature evolution around the one-phonon threshold, where Polaron Interband Transitions appear. 



u> w 2Ep. Clearly, the transfer of spectral weight seen 
in Fig0]a when going from T = 0.03w to T — 0.2w 
must be traced back to the change of behavior in the 
excitations at v < Eo described here. 

A more careful look at Fig. [S]a (expanded in Fig. 
EJb) shows that the change in the overall distribution 
of ImEe _ ^" _£;o ' is accompanied by a slight shift of the 
position of the individual peaks. In fact, as can be shown 
by direct inspection of the continued fraction expansion 
of Ref. 7 ., each peak at v < Eq is a replica of the polaron 
band at Eq, of width W . At temperatures T <C W, due 
to the exponential weighting factor, only the low mo- 
mentum states close to the bottom edges contribute (cf. 
Fig. Ob). Such states are free-electron like, and give 
rise to an optical absorption spectrum with a weak cou- 
pling character. As T > W, on the other hand, more 
localized higher momentum states come into play, and 
the polaronic behavior is recovered. We remark that this 
"reentrant" behaviour is not restricted to the zero den- 
sity case treated here. It can be found also in a two site 
cluster— as well as in the Holstein model at half filling 30 . 
The reason is ultimately due to the different roles played 
by the quantum and thermal fluctuations in the polaron 
crossover region. Quantum fluctuations are known to 
stabilize the non-polaronic phase, while incoherent fluc- 
tuations such as static disorder or thermal fluctuations 
stabilize the polaron. As the temperature increases, po- 
larons are therefore first stabilized, before they eventually 
dissociate at a higher temperature. 

According to the above arguments, the transfer of spec- 
tral weight illustrated in Fig. 0Ja is ultimately controlled 
by a temperature scale set by the renormalized band- 
width W, whose evolution with the interaction strength 
is shown in Fig. In the strong coupling regime, where 
W vanishes exponentially, such phenomenon can hardly 
be observed in practice. On the other hand, the very ex- 
istence of a polaron subband separated from the contin- 



uum of excited states requires moderately large values of 
A, and the coexistence of electron-like states and localized 
states is specific to the adiabatic regime. The competi- 
tion between all these conditions explains why the low- 
temperature transfer of spectral weight described here 
occurs in the vicinity of the adiabatic polaron crossover, 
roughly for A c — O.57 < A < A c + 7 and 7 < 0.5 (shaded 
area in Fig [21 . 

A formal expression for the optical conductivity in the 
limit T <C W is derived in App. B, that we reproduce 
here: 

cr(w) = [1-W(0)]- x (14) 
EjU / deN(e)<j ) (e)p> (e, lj + E Q — puj Q )u p {e) 

where W(0) is a constant between and 1 depending 
on the interaction parameters, and the functions u p (e) 
are defined in Eq. I|62|) . This formula replaces Eq. (7) 
of Ref. 18 , which was incorrectly identified with the zero 
temperature limit of the optical absorption but is instead 
a preasymptotic contribution as T — > 0. 

B. Polaron Interband Transitions 

A detailed view of the low frequency optical absorp- 
tion in the intermediate coupling regime is presented in 
Fig. 0Jb. As was mentioned in the previous Section, 
the spectrum at the lowest temperature (full red curve) 
has a sharp threshold at w = luq, followed by an ab- 
sorption pattern which roughly reproduces the shape of 
the single particle spectral density, shifted by loo [see the 
weak coupling formula (|25|l in the appendix, as well as 
its generalization Eq. (|14|) ]. For example, the peak at 
u> ~ 0.17 is directly related to the maximum observed 
in the spectral density at the top of the polaron band 
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FIG. 5: (Color online) a) Spectral density iV* and weighted 
scattering rate ImY,(v)e~ ( v ~ Bo) at A = 0.9 and 7 = 0.1, at 
different temperatures, b) Expanded view of a generic n-th 
peak in the weighted scattering rate (n = 3, indicated by an 
arrow in the main panel, whose width equals the renormal- 
ized polaron bandwidth W (the data at T = 0.03u>o have been 
multiplied by a factor 10 3 ). The weight moves from the bot- 
tom to the top of the band at T ~ W . c) Expanded view of 
the polaron band in the spectral density N* , separated from 
the high-energy continuum by a gap. 
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FIG. 6: (Color online) Width of the polaron band as a func- 
tion of the interaction parameter A for different values of the 
adiabaticity parameter 7, at T = 0. 



(compare Fig. 0]b with Fig. 0c), and its position co- 
incides with the expected value to = uj + W — 0.173. 
Upon increasing the temperature, the most striking effect 
in this low-frequency region is the emergence of a sharp 
asymmetric peak within the gap, whose intensity follows a 
puzzling non-monotonic temperature dependence. Such 
resonance has already been reported in Ref^&, where it 
was termed polaron interband transition (PIT) , and sug- 
gested as a possible explanation of the peaks experimen- 
tally observed in Refi^i. In order to understand its ori- 
gin, one has to go beyond the low-temperature arguments 
presented in the previous Section, and take into account 
the contributions to the convolution integral in Eq. (JSJ 
coming from the states at v > Eq. 

Because of thermal activation, there is a non- vanishing 
probability for transitions involving an initial state at en- 
ergies above the ground state. The Fermi golden rule tells 
us that the transition probability at a given frequency uj 
will be proportional to the conjugate density of the ini- 
tial and final states separated by uj, times the occupation 
factor of the initial state. Therefore, a sharp resonance 
will arise if: i) many pairs of initial and final states exist 
that are separated by approximately the same energy uj* , 
which is equivalent to the condition that the energy dis- 
persions in the initial and final subband are parallel; ii) 
such states are long lived, which is possible at low tem- 
peratures because JraS is exponentially suppressed in 
the whole region Eq < v < Eq + luq . Both conditions can 
be fulfilled in the intermediate coupling regime, as illus- 
trated in Fig. 3 of Ref^. Upon increasing the tempera- 
ture, the states involved in the transitions are eventually 
smeared out by thermal disorder, which explains the non 
monotonic temperature dependence of this resonance. 35 

This qualitative argument can be formalized as follows 
(the detailed calculations are presented in Appendix B). 
When the scattering rate in the initial state is negligible, 
the spectral function p(e, v) can be replaced by a delta 
function and the convolution integral in Eq. (jSJ) becomes 
[cf. Eq. gT}]: 



V(uj) 



dvN{v - ReL{v))<j){v - ReL{v))e- p{ - v - Eo> > 



B 



x Im 



1 



7T uj + ReY,(v) — ReY>{v + lu) +iT 



.(15) 



This quantity is maximum for those values of uj such that 
there is a wide range of energies v where the denominator 
is close to zero. This condition corresponds to the pair 
of equations 



ReYl{v) - ReY;(u 



dReY* 



dReY, 



Oh 



- uj) - 
= 







(16) 



whose solution we denote as uj* and v* . The optical 
conductivity can be evaluated by expanding Eq. fl"5|l 
around these special values, and taking T = — JmS(i/ + 
uj) to be small and constant in the frequency range of 
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interest, which gives 

TT (1 - e-P")e-P v ' 



(jpit{w) 



(17) 



xN(u* - ReE(v*))</>(v* - ReT,{v*))T 



M(u-u*) 
f 



where M is the effective mass of the initial (and final) 
states and Q is related to the curvature of the dispersion. 
The function 



Hv) = 



\jy + ^y 2 + ^ 



(18) 



has an inverse square root tail ~ y^ 1 ^ 2 at large y, and 
tends to 1 for < y < 1, giving rise to a narrow peak 
of width Alo ~ T/M. Such asymmetric shape can be 
clearly identified in FigQ]b, at uj* ~ 0.3wo- Reminding 
that r ~ e~^ w °, and using the definition Eq. @, we see 
that the weight of the narrow peak scales as 



e -l3uj* £ -[3(v*-E +uo/2) 



(19) 



This quantity rises exponentially with temperature and 
then decreases as a power law at temperatures above 
some fraction of whose precise value depends on the 
parameters u>*, v* and W. 



VI. CONCLUSION 

In this work, we have studied the optical absorption 
of electrons interacting with phonons in the framework 
of the Holstein model, applying a unified treatment - 
the dynamical mean-field theory — which is able to ac- 
count for both the quantum nature of the phonons and 
the effects of a finite electronic bandwidth. The present 
non-perturbative approach allows to span all the ranges 
of parameters of the model, from weak to strong electron- 
phonon interactions, and from the adiabatic to the an- 
tiadiabatic limit. The limits of validity of the standard 
available formulas are pointed out, and the quantitative 
and qualitative deviations arising in wide ranges of the 
parameters space are analyzed, which will hopefully lead 
to a better understanding of the experiments performed 
on polaronic systems. 

In the antiadiabatic regime, the present method al- 
lows to address the evolution of the optical absorption 
away from the single molecule limit, where the spectra 
reduce to a distribution of delta functions. Consider- 
ing a finite bandwidth D yields an intrinsic width to the 
multi-phonon peaks, which is not uniform over the whole 
spectral range (peaks at high frequency are broader than 
at low frequency) and generally does not obey the expo- 
nential narrowing predicted by the standard approaches. 
Moreover, in the strong coupling regime, the peaks are 
found to broaden upon increasing the temperature, in 



marked contradiction to the commonly accepted results 
in the literature, which are based on uncontrolled approx- 
imations. 

The situation is different in the adiabatic regime, where 
it is found that the nature of the optical absorption 
changes drastically at the sharp polaron crossover oc- 
curring at A ~ A c . The weak coupling formula i|l(J|) is 
shown to extend qualitatively to finite values of A, up 
to the polaron crossover, whose proximity is signaled by 
the emergence of a fine structure related to multi-phonon 
processes. Beyond the polaron crossover, the absorption 
mechanism in the strong coupling regime is found to de- 
pend on an additional parameter, the ratio s/ D between 
the broadening of the Franck-Condon line and the free- 
electron bandwidth. The usual polaronic gaussian line- 
shape described by Eq. (12) is recovered for s ^> D. In 
the opposite limit, s <C D, a different absorption mech- 
anism sets in, related to the photoionization of the po- 
laron towards the free-electron continuum. Correspond- 
ingly, the lineshape in this limit is entirely determined 
by the shape of the free-electron band, and is charac- 
terized by a sharp absorption edge at finite frequency, 
as described by Eq. (13). The optical spectra at inter- 
mediate values of s/D, a situation often encountered in 
real systems, are not properly described by any of the 
two limiting formulas Eqs. (12) and (13). In particular, 
due to finite bandwidth effects, the frequency of the ab- 
sorption maximum is appreciably lower than the usual 
estimate uj max — 2Ep, and the lineshape is more asym- 
metric than the gaussian of Eq. (12), which should be 
taken into account when interpreting experimental data. 

In the intermediate coupling region around A ~ A c , 
qualitatively new features arise that are distinctive of 
the adiabatic polaron crossover. First, the optical ab- 
sorption exhibits a reentrant behavior, switching from 
weak-coupling like to polaronic-like upon increasing the 
temperature. The temperature scale that governs such 
transfer of spectral weight is of electronic origin, being 
set by the renormalized polaronic bandwidth W, a quan- 
tity that is necessarily less (or equal) than the phonon 
energy, and rapidly decreases upon increasing A. In ad- 
dition, sharp peaks with a non-monotonic temperature 
dependence emerge at typical phonon energies, and are 
most clearly visible in the region below the single-phonon 
absorption gap. These can be ascribed to thermally ex- 
cited resonant transitions between long lived states lo- 
cated in different subbands in the polaron internal struc- 
ture, and their peculiar temperature dependence results 
from the competition between the thermal activation and 
the thermal broadening of the corresponding states. 

To conclude, let us observe that the present formalism 
can in principle be applied to any finite-dimensional lat- 
tice by replacing N(e) of Eq. (7) with the appropriate 
DOS. It has already been mentioned (cf. end of Section 
V A) that Eq. (13) reproduces very accurately the one- 
dimensional exact diagonalization data of Ref. 13 obtained 
for s <C D. Even away from this limit, however, prelim- 
inary DMFT results obtained using the one-dimensional 
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DOS N(e) — l/(ir\/D 2 — e 2 ) are in qualitative agree- 
ment with the numerical results of Refi^, at least in the 
high frequency part of the spectra. In particular, in the 
intermediate coupling regime, the polaronic absorption 
maximum appears to be shifted downwards as compared 
with the three-dimensional case (cf. Figs. 4 and 6 in 
RefiiS). This result can be understood from the argu- 
ments presented in Section V A by observing that due to 
the singularity of the one-dimensional DOS, the function 
N* appearing in Fig. is peaked around v — Eq ~ 
instead of v — E — Ep. 

Concerning vertex corrections, which are absent in our 
formalism, we have seen that these are irrelevant in the 
limiting regimes of weak and strong coupling, as well as in 
the adiabatic regime s -C D at any coupling strength and 
for any lattice dimensionality (including in one dimen- 
sion). On the other hand, their presence could quantita- 
tively change the optical spectra when both the phonon 
frequency and coupling strength lie in the intermedi- 
ate regime. From this point of view, the polaron inter- 
band transitions evidenced in the intermediate adiabatic 
crossover region, relying on the resonance condition Eq. 
(|16fl . could indeed be affected by the inclusion of vertex 
corrections. 

As a last remark, we stress that the description of the 
optical properties provided in this work, based on the so- 
lution of the Holstein model for a single particle, is valid 
in principle for a system of independent polarons. How 
the above features are modified by polaron-polaron in- 
teractions at finite densities remains a fundamental ques- 
tion, that needs to be addressed in the future. 



Appendix A: Derivation of limiting formulas 
Weak coupling limit, low temperatures T -C wo- 
rn the limit of weak interactions, we can write the spec- 
tral function p e appearing in Eq. JSJ as 



Pe = P. 



(o) 



Spe 



where pi°^ is non-interacting and reduces to a delta func- 
tion. Substituting the previous expression in Eq. (J5J and 
taking into account only the contributions proportional 
to Sp we get 



2?(w) = / deN(e)<P(e) J dve-^ v ~ E ^ x (20) 
x pf){v)5p e {p +oj) + 8p e (v)pQ\v + u) 
Performing the two integrals in v yields 
V(u) = \ fM(e)#)[e-*- E » ) /raE(e + i l i)+ 



JmE(e — u>)e 



-(3(e-u>-E ) 



Perturbation theory gives, to order g 2 and at zero density 

E(w) - g 2 [(1 - n B )G{u - cj ) + n B G(u> + cj )\ (22) 

where n B is the Bose occupation number. Substituting 
the self-energy into Eq. l(2*T|l gives rise to 4 terms, namely: 

V{oj) = ^2 / deN(e)4>(e) x 

x (l-n B )e- p(e ~ Eo) N(e + uj-uj ) (23a) 

+ n B e-P {e ~ Eo) N(e + u} + u} ) (23b) 
+ (1 - n B )e- /3< - e - u '- Eo) N(e - to - oj ) (23c) 



n B e 



-0(e- 



(23d) 



At temperatures T <C u>o, the second and third term 
are exponentially small, and can be neglected. Note 
also that, since the above terms are explicitely propor- 
tional to g 2 , the normalization factor JBJ can be evalu- 
ated in the non- interacting limit, which gives Af((3) = 
{2/f3)e- f3D h(l3D), where h is the modified Bessel func- 
tion of the first kind, whose limiting behaviors are J\f ~ 
y/2/^/(PD)^ 2 for (3D > 1 and M -> 1 for (3D < 1. 

At temperatures much smaller than the width of the 
band under study, due to the presence of the exponential 
factors, only the low energy edges of the density of states 
(corresponding to low-momentum states) contribute to 
the integrals in Eq. In this case we can make use 

of the following approximation: 

dee-^\e-vff{e) ~ ^T(t + 1) (24) 

where T(n) is the Euler gamma function and f{y) is a 
generic continuous function. By substituting the previ- 
ous expression in equations only the term (|2*3T 
survives in the limit T -C D, and we obtain 



g 2 ir 



= ^</>(w -D- lo )N(uj -D- uj ). (25) 



which corresponds to a sharp edge behavior (t(uj) ~ 
(luq — w) 3 / 2 , in agreement with perturbative calculations 
in three dimensions^. Similar arguments lead to <j(uj) ~ 
{uj -uj) d / 2 if one assumes the appropriate d— dimensional 
noninteracting band in Eq. (J5J . 

In the adiabatic regime, the condition T <C wq ensures 
that T <C D. In the opposite antiadiabatic regime, the 
result remains valid for T < D, while for D < T < 
u>o, we can set the exponential factors equal to 1 in Eqs. 
(|?3l . leading to 



(26) 



where the function F is defined as 



(21) 



r-D-lvl 

F(y) = I dxN(x)N{x- 



\y\)[<t>{x)+c}>{x+\y\)] (27) 
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(note that current vertex corrections are not included 
in this formula; however, they lead to the same func- 
tional form in both the adiabatic and anti-adiabatic lim- 
its). In this temperature range, the absorption related to 
the excitation of a single phonon is nonzero in the range 
luq — 2D < lu < lu + 2D and is symmetric around lu = luq 
(see the inset in Fig. Q]a) where a = (128/45ir)g 2 D/uj$. 
By comparing Eq. (|27|) with Eq. I|25|) , we see that in the 
antiadiabatic regime, the absorption peak is shifted to 
lower frequency when the temperature is increased above 
T ~ D, and the sharp square root edge is converted into 
u{w) ~ [w - (lu - 2£>)] 3 . 



Strong coupling limit, narrow bands. 

In the adiabatic limit (ojq — > 0), the strong coupling 
regime can be attained by increasing A while keeping 
the polaron energy Ep finite. In this case, the spec- 
tral function of a single polaron becomes independent of 
momentum (the e index drops because the bandwidth 
D = Ep/X — > 0). The self-energy in Eq. (2) becomes 
£ = —gX, where X is a classical variable that acts as a 
source of disorder that modifies the electronic levels ac- 
cording to the Boltzmann distribution^^. Correspond- 
ingly, the spectral function tends to a gaussian 



p(u) = 



1 



V27TS 2 



exp 



v 

2J 2 



(28) 



whose variance s = \J2EpT is determined by the thermal 
fluctuations of the classical phonon field. Once the e- 
integral is factored out (yielding a numerical prefactor), 
the convolution integral in Eq. JSJ is readily evaluated 
and yields 



cr(w) 



4lu y/Ai 



exp 



(lu - 2E P ) 

4^ 



(29) 



which coincides with the standard result of Refs. 5 ' 9 ' 10,11 
It should be noted that, due to the initial assumption 
lu$ = 0, the above derivation is strictly valid only in the 
high temperature limit T > wo- However, Eq. (|29(l can 
be generalized to all temperatures by setting 



s 2 (T) = E p lu coth(^) 



(30) 



[note that Eq. I|28|) is also valid at T = 0, in which 
case the variance is determined by the phonon zero point 
fluctuations and becomes s — V TSpujgLZS^. Eq. is 
valid in the strong coupling adiabatic regime provided 
that s > D. 

Note that in the limit of an isolated molecule (D = 0) 
at finite luq, which is qualitatively representative of the 
antiadiabatic regime, the spectral function consists of 
a distribution of delta functions, describing the multi- 
phonon excitations inside the polaron potential- well 5 . 



The convolution integral in Eq. (|SJ) can be carried out, 
and yields a discrete absorption spectrum, whose enve- 
lope in the strong coupling limit a 2 = Ep/luq S> 1 coin- 
cides with Eq. l|2^l) . 



Strong coupling limit, wide bands. 

There is an alternative way of reaching the polaronic 
adiabatic regime, i.e. assuming that the phonon induced 
broadening of electronic levels s is negligible compared 
to the electron dispersion D. For sufficiently large Ep, 
the optical absorption will be due to transitions from a 
perfectly localized level whose electronic energy is —2Ep, 
to the free-electron continuum. Put in the context of 
the convolution integral Eq. JSJ, this corresponds to the 
following replacements 

p(e,Lu + v) = 8(e — lu — v) (31) 
p{e,v) = 5{v + 2E P ) (32) 

where the last quantity is independent of e. The result is 



cMs«-D«-Ep = 7T- 



(lu-2E p )N{uj~2Ep) (33) 



This result is equivalent to that reported in Refii2ii^ 
for small polarons in one-dimensional systems, and in 
Ref^i for large polarons in two and three dimensions. 
However, the limit of a perfectly localized state assumed 
here is strictly valid only for Ep — > oo. At finite Ep, a 
more accurate formula can be obtained by restoring the 
e-dependence of the initial state in Ea. (|32|l . This can be 
achieved by taking 



1 1 1 imEfr/) 
P(e, v) = Im 7 ^ 



(34) 



where we have used the fact that both ReT, and JmS are 
small around v ~ —2Ep. The scattering rate at these 
energies can in principle be extracted from the atomic 
limit of Refi 5 -, because it is a local quantity. For s — * it 
is given by: 

- -ImY,{v)e- f3{v+Ep) = v 2 5{v + 2E P ) (35) 

Evaluation of Eq. (JjjJ now yields a more asymmetric 
shape 

a(w) s « D = 7r 4£ ; pl ~ e - - 2E p )N(lu - 2E P ) 

UJ Z LU 

(36) 

which reduces to the previous Eq. JsSJ in the limit E P > 
D. 

It should be noted that when the appropriate one- 
dimensional DOS is used, Eq. (|36|l describes much better 
the exact diagonalization data of Ref. 13 than Eq. (|33l) it- 
self, that was used by Alexandrov and coworkers in their 
Figs. 3.cand3.d. A possible generalization of the present 
scheme to finite values of s is currently under study^ 
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Appendix B: Optical conductivity at low 
temperature 

In this appendix we give a formal expression for the op- 
tical conductivity in the limit of low temperatures, which 
is useful for the understanding of the results in the inter- 
mediate coupling regime. 

We can formally separate p(e, v) in two parts, at ener- 
gies above and below the ground state 

p(e, v) = 6{v - E )p> (e, v) + 9(E - u)p< (e, v) (37) 

[a similar separation will hold for N*{v)] Accordingly, 
in Eq. Q we can separate in both numerator T> and 
denominator Af the contributions coming from energies 
lower and higher than the ground state energy E , ob- 
taining 



V >> (lu)+V <> (lu)+V << {lu) (38) 
Af > +Af < (39) 



In Eq. (|38|) the last contribution is negligible at low tem- 
perature for lu > since p < (lu + v) decays exponentially 
as exp[— /3(Eq + ui — v)]. On the contrary both Af < and 
V <> (lu) cannot be neglected due to the thermal popula- 
tion factor exp[— fi(v — Eq)] appearing in Eq. Q. Let us 
calculate the remaining four terms separately. 

Calculation ofT> >> (lu) 

/>oo pD 

V»{u) = dv deN(e)(f>(€)e-^ u - Eo) x 

J Eg J-D 

x p>(e,v)p>(e,LU + v) (40) 

At low temperature, p > (e, v) weighted by the exponential 
factor gets contributions only for v ~ Eq. This occurs 
when the temperature is so low that both the thermal 
occupations of the first excited subband and of the inco- 
herent background edge at v = Eq + ujq are negligible, 
i.e. for (3luq ^> 1. In this case the scattering rate is neg- 
ligible (the low-energy states are coherent), p > (e 1 v) can 
be replaced by 5[v — e — ReYj(v)\ and the integral in e 
appearing in Eq. (|40|) can be carried out: 



D >> (w) = / dvN{v - ReE{v))4>{v - ReE(v)) 

Je 

x e~ l3( ' J - Eo) p > {v-ReY,{v) 1 La + v). (41) 

At temperatures lower than the renormalized bandwidth 
W, the thermal weighting factor selects an integration 
range of width T around the ground state energy where, 
from Eqs. J7|) and ©, we have 



N(u-BeS(u)) ~ -^^2D/Z{ v -EqY/ 2 (42) 



(y-ReS(y)) ~ ^J^-(v - E ) 



(43) 



Here Z = (l — d^ReY.) 1 is the ground state quasiparticle 
residue, equal to the inverse of the effective mass in the 



case of a purely local self-energy. Making use of Eq. I|24() 
and the definition Eq — ReT,(Eo) = —D, we obtain 



V»(w) = C(T) 



—Im — 

7T LU 



Eq+D-Z(lU + Eq) 



_ (44) 
with C(T) = D 2 /Z 3 / 2 [T/Df/ 2 for T < W. In 

the strong coupling regime, the polaron band becomes 
extremely narrow and can itself be replaced by a delta 
function. For T 3> W — > 0, we obtain the same expres- 
sion Eq. (glj with C(T) ~ iW 2 — > 0, and this contribu- 
tion becomes negligible. 

Polaron interband transitions The previous formula 
for T> >> was derived assuming that T < wq. Upon in- 
creasing the temperature, other contributions arise that 
involve states at frequencies v away from the ground state 
Eq . These can give rise to sharp resonances in the optical 
absorption spectra (see Section IV B). Eq. (|41|1 can be 
rewritten as 

/>oo 

■D»(w)= dvN{v-Re£{v))4>(v-Re?,{v)) (45) 



xe 



-0(v-E Q ) 



— I 111 

7T (j + ReS(u) - ReE{v + u) + iT 



where T = —ImYi(v-\-Lu) is taken to be small and constant 
in the frequency range of interest. A resonance will arise 
for those values of lu such that there is a wide range of 
energies v where the denominator is close to zero. This 
condition corresponds to the pair of equations 



ReT,(v) - ReY,{v + lu) = 



dReY, 



dv 



dReY, 



dv 







V+UJ 



whose solutions we denote lu* and v* . Now the denomi- 
nator in Eq. (|4B|) can be expanded around these special 
points, which gives 



M{lu 



iT 



M = 1 — (dYi/dv) u * being the (equal) effective mass of 
the initial (and final) states and Q a parameter related to 
the curvature of _Re£. If one assumes that the prefactors 
N and 4> are smooth functions in the vicinity of v* , these 
can be taken out of the integral, which can be evaluated 
to 



V 



PIT 



N{v* - ReH{v*))4>{v* - ReZ{v*)) 



xe 

where the function 



-^*jF[M(w-w*)/r] 



Vv 2 



i 



(46) 



(47) 



tends to 1 for < y <C 1 and has an inverse square root 
tail ~ y~ x l 2 at large y. 
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Calculation o/2? <> (tj) 
V<>(u) = I dv I deN(e)cf>{e)e~^~ E ^ x 

JEo-LU J-D 

x p<(e,v)p>(e,u + u) (48) 

The excitations in p < are mostly incoherent, being due 
to the presence of thermally activated phonons. This is 
best expressed by showing explicitely the proportionality 
to the scattering rate 



1 



P < (e,v)= — 



ImY>(v) 



n(u- fte£(» - e) 2 + (JmE(v)) : 



(49) 



Through Eq. jI5|L the product e^'^^'Mfi/) ap- 
pears in Eq. 1(48(1 . A sample plot of this product is re- 
ported in Fig. 0b. At energies below the ground state, 
taking advantage of the continued fraction expansion of 
Refi, the imaginary part of £ can be approximately writ- 
ten as a sum of separate contributions 



P =i 



u p {v) 



(50) 



As was the case for T> >> , we treat separately the cases 
T <C VF and T 3> W — > 0. In the former case, u p is a 
function with a square root edge at v — Eq — pujQ 

u p (y) = u p (y - E a + pw ) 1/2 



(51) 



which is reminiscent of the behaviour of the spectral den- 
sity at energy Eq — (p— 1)o;q and u p is an unknown weight- 
ing factor. As can be seen from Fig. 0b, each u p {v) has 
a bandwidth which is independent on p, and equal to the 
width W of the first polaronic band. Therefore, when 
T <C W, the square root edges can all be replaced by 
delta functions by making use of Eq. 1(241) . When Eqs. 
(|50I51|I are substituted in Eq. H48|) the z/-integral can be 
carried out, leading to 



v <:> (lu) = r(3/2)r 



:;/2 E 
P =i 



D 



deiV(e)0(e) x 



x p>(e,u + E -pu! ) x 

1 Un 



7r {E Q - pu - e - ReT,(E - pw )) 



where p — int(uj / uj ) . In Eq. itB^ we have neglected 
the exponentially vanishing scattering rate JmE in the 
denominator of Eq. (|49l) . At higher temperatures, or in 
the strong coupling regime, where ujq T » W — > 0, 
we can take 



u p (v) = Up5(^ - E Q - W + ploq) 



(53) 



which leads to an expression analogous to Eq. 1)52(1 , with 
the prefactor replaced by e~@ w and Eq — > .Eq + VF, re- 
flecting the fact that the dominant weight is now carried 
by localized excitations at the top of the band (see the 
discussion in Sec. Ill A). 



Calculation of Af' 



Af> = 



dv 



deN{e)e- p[u - Eo) p>(e, V ). (54) 



Since a gap of width loq separates the coherent pole of 
p > from the incoherent continuum, for (3ujq 3> 1 only the 
lowest energy edge of the spectral density contributes to 
the integral Eq. Using Eqs. 124142(1 we get 



Af' 



ttD 2 



y/2D/ZT(3/2)T 



3/2 



(55) 



Calculation of Af < 
Af< -- 



E pD 

dv j deN(e)p < (e,v). (56) 

•oo J-D 



As was stated before, p < represents incoherent states as 
given by Eq. (|50|) . For f3W 1, the integrals in v can 
be carried out with the help of Eq. 1(24(1 giving 



OO .£) 

r(3/2)T 3/2 V / deN{e) 
P =i J - D 



Af < 



7r (E - pu - e- ReY,(E Q - _pw )) 2 



(57) 



At higher temperatures, T » W, we get an analogous 
result, with the prefactor replaced by e~@ w and Eq — > 
Eq + W. By comparing the two expressions, it can be 
deduced that the term Af K changes behavior from ~ T 3 / 2 
to e'PW at a temperature T w ~ W/log(D/W). The 
resulting normalization factor Af has the form ~ T 3 / 2 in 
the weak coupling regime and ~ e~@ w — > 1 in the strong 
coupling regime. In the intermediate coupling regime, 
where W is of the order of u)q, a crossover between the 
two behaviors occurs at Tw- 

Final formula for o~{ui) at T — 0. Collecting the re- 
sults Eqs. ((44I52I55I57|) into Eq. Q) we get the final 
result. To separate the terms (|44f) and 1(55(1 which in- 
volve the coherent states close to the ground state Eq 
from the remainder, it is useful to introduce a tempera- 
ture dependent weight 



W{T) 



Af> 
Af 



(58) 



which tends to a constant for T <C W, and to a (different) 
constant for T > T w . Then 

a(uj) = W(T)cro(w) + [1 - W{T)} <ti(u>) (59) 

where, from Eqs. ((44I55|) : 

1 \ 



. . ttD/Z ( 1 
a (uj) = — T I Im 



7T uj + Eq + D - T,(Eq + uf) 



(60) 
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Eq. H60(l is Eq. (7) of Ref. 18 , where it was incorrectly 
identified with the zero temperature limit of the optical 
absorption (note that the linear dependence on temper- 
ature makes this contribution vanish at T — 0). 

Under the above hypothesis the remaining term o\ 
does not depend explicitely on temperature. In fact, us- 
ing Eqs. i|52l57[l we get a formal expression 

_ 7T El=i I deN{e)<f>(e)p> (e, u + Eq - P io )u p (e) 
- E™ildeN(e)u p (e) 

(61) 
(62) 



where 

u p (e) 



(E - puj - e - ReE(Eo - puo)) 2 



and at T = we have a(oj) = [l-VV(0)]oi(w). Note that, 
since the sum is limited to p, at low temperature only the 
first term in the numerator of Eq. I|61|l contributes to the 



Drude peak at oj ~ and to the one-phonon threshold 
at u> = uj . 

To recover the perturbative expression Eq. I|25|l we 
consider only a single contribution in the sums appear- 
ing in Eq. I|61|l . Then we perform the integral in the 
numerator of Eq. (|61|) using the free spectral function 
p > (e, lu + Eq — pu>o) — S(lo + Eq — pujQ — e) obtaining 



cri(w) 



7T N(ld + Eq — UJq)4>{lO + Eq - UJq)Ui(u + E — OJq) 



To proceed further we have from Eq. 122(1 

,2\/2~D 



ui(a> + E Q - cjq) 



D 2 uj 2 



(63) 



(64) 



where we have neglected the self-energy term in the de- 
nominator of Mi for (3ojq ^S> 1. Eq. (|25|1 follows by multi- 
plying (ESI) by the factor (1 - W). 



1 I.G. Austin and N. F. Mott, Adv. Phys. 18, 41 (1969) 

2 T. Holstein Ann.Phys. 8, 325 (1959), ibid. 343 (1959). 

3 S. Fratini and S. Ciuchi, Phys. Rev. Lett. 91, 256403 (2003) 

4 A.S. Alexandrov and J. Ranninger, Phys. Rev. B 45, 
R13109 (1992); J. Ranninger Phys. Rev. B 48, R13166 
(1993) 

5 G.D. Mahan, Many-Particle Physics, Plenum Press, New 
York, 2nd edition (1990) 

6 H. Sumi, J. Phys. Soc. Jpn. 36, 770 (1974) 

7 S. Ciuchi, F. de Pasquale, S. Fratini and D. Feinberg, Phys. 
Rev. B 56, 4494 (1997) 

8 V. Perebeinos and P.B. Allen, Phys. Rev. Lett. 85, 5178 
(2000) 

9 D. M. Eagles, Phys. Rev. 130, 1381 (1963); M. I. Klinger, 
Phys. Letters 7, 102 (1963); V. N. Bogomolov, E. K. Kudi- 
nov, D. N. Mirlin and Yu. A. Firsov, Fiz. Tverd. Tela 9, 
2077 (1967) [Sov. Phys. Solid State 9, 1630 (1968)] 

10 H. G. Reik. Sol. St. Comm 1, 102 (1963); H. G. Reik and 
D. Heese, J. Phys. Chem. Solids 28, 581 (1967); H. G. Reik 
in " Polarons in Ionic Crystals and Polar Semiconductors" , 
ed. J.Devreese, North-Holland, Amsterdam (1972) 

11 D. Emin, Phys. Rev. B 48, 13691 (1993) 

12 Yu. A. Firsov, Fiz. Tverd. Tela 10, 1950 (1968) [Sov. Phys. 
Solid State 10, 1537 (1969)] 

13 A. S. Alexandrov, V. V. Kabanov and D. K. Ray, Physica 
C 224, 247 (1994) 

14 I. G. Lang and Yu. A. Firsov, Zh. Eksp. Teor. Fiz. 43, 1843 
(1962) [Sov. Phys. JETP 16, 1301 (1963)] 

15 V. V. Kabanov and O. Yu. Mashtakov, Phys. Rev. B 47, 
6060 (1993). 

16 M. Capone, S. Ciuchi, and C. Grimaldi Europhys. Lett. 
42, 523 (1998) 

17 G. Wellein and H. Fehske, Phys. Rev. B 58, 6208 (1998) 

18 S. Fratini, F. de Pasquale, and S. Ciuchi, Phys. Rev. B 63, 
153101 (2001) 

19 G. Schubert, G. Wellein, A. Weisse, A. Alvermann, and H. 
Fehske, Phys. Rev. B. 72, 104304 (2005) 

20 A.Georges, et al, Rev. Mod. Phys. 68, 13 (1996) and ref- 
erences cited therein. 

21 E. Jeckelmann, H. Fehske, |cond-ma t/0510637 



22 A. Khurana Phys. Rev. Lett. 64, 1990 (1990). 

23 W. Chung and J.K. Freericks, Phys. Rev. B 57, 11955 
(1998) 

24 A. Chattopadhyay, A.J. Millis, S.Das Sarma, Phys. Rev. 
B 61, 10738 (2000) 

25 N. Bliimer and P. G. J. van Dongen, in "Concepts in Elec- 
tron Correlation", Eds.: Alex C. Hewson, Veljko Zlatic, 
NATO Science Series, ISBN 1-4020-1419-8, Kluwer (2003) 
|cond-mat /0303204 

20 S. Fratini and S. Ciuchi. Phys. Rev. B 72. 235107 (2005) 

27 M. Capone, W. Stephan, M. Grilli, Phys. Rev. B 56, 4484 
(1997) 

28 M. J. Rozenberg, G. Kotliar, H. Kajueter, Phys. Rev. B 
54, 8452 (1996) 

29 S. Paganelli, S. Ciuchi, cond-mat/0410032 

30 S. Ciuchi, G. Sangiovanm and M. Capone, 
cond-mat/0503681 

31 Y. H. Kim, A. J. Heeger, L. Acedo, G. Stucky, F. Wudl, 
Phys. Rev. B 36, R7252 (1987); A.V. Bazhenov et al, 
Physica C 214, 45 (1993); P. Calvani et al., Europhys. Lett. 
31, 473 (1995); P. Giura, Ph.D. thesis, Universite Paris XI 
Orsay (2000) (unpublished) 

32 S. Fratini, S. Ciuchi, in preparation 

33 It is known from the atomic limit that, once multiplied by 
e -P(v-E ) ^ S p ec t ra i W eight at negative frequencies has 
a distribution centered at v — Eo ~ — Ep — see Refi-, and 
Fig. El 

34 Additional peaks have also been reported to arise in 
the adiabatic intermediate coupling regime, in one space 
dimension 19 ". Such peaks have been ascribed to transitions 
between polaron states with different spatial structures, 
and cannot be described by the present (local) DMFT 
treatment. 

35 More generally, when more than one subband is present 
in the excitation spectrum, a fraction or even all of the 
corresponding states can have a vanishing scattering rate 
at low temperatures (cf. Figs. 11 and 13 of Reft->), and 
several PITs can appear. 



